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The study of optomechanical systems has attracted much attention, most of which are concentrated in the physics in the small- 
amplitude regime. While in this article, we focus on optomechanics in the extremely-large-amplitude regime and consider both 
classical and quantum dynamics. Firstly, we study classical dynamics in a membrane-in-the-middle optomechanical system in 
which a partially reflecting and flexible membrane is suspended inside an optical cavity. We show that the membrane can present 
self-sustained oscillations with limit cycles in the shape of sawtooth-edged ellipses and exhibit dynamical multi stability. Then, we 
study the dynamics of the quantum fluctuations around the classical orbits. By using the logarithmic negativity, we calculate the 
evolution of the quantum entanglement between the optical cavity mode and the membrane during the mechanical oscillation. We 
show that there is some synchronism between the classical dynamical process and the evolution of the quantum entanglement. 
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1 INTRODUCTION 

Optomechanics m concerns with the coupling between op¬ 
tical cavity modes and mechanical degrees of freedom via 
radiation pressure (3j, optical gradient forces 0. photother- 
mal forces [5] or the Doppler effect f6|. Optomechanical sys¬ 
tems have potential applications in many areas, such as detec¬ 
tion of gravitational wave 00’ sensitive sensors ]T0l}T2) , 
mechanical memory (T3), cooling of mechanical resonators 
G3G3. coherent manipulation of light |]T6}[T7), studies of 
quantum entanglement | T8J-2TJ, preparation of macroscopic 
quantum state 122 ]24| and quantum information processing 


Most early works focus attention on optomechanics in the 
small-amplitude regime, in which the optomechanical cou- 
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pling can be approximately considered as linear or quadratic. 
With a driving laser of frequency blue detuned with respect to 
the cavity resonance and power above a certain threshold, the 
mechanical resonator can run into self-sustained oscillations 
f27l(29] | and dynamical multistability may emerge |[30]|3 1 [ ■ 
The mechanical resonator conducts approximately sinusoidal 
oscillations, and therefore, its limit cycles in the phase space 
are approximately elliptical. 

In the large-amplitude regime, which can be realized in an 
optomechanical system driven by a high-power laser, the me¬ 
chanical resonator can display different self-sustained oscil¬ 
lations with limit cycles being mushroom-like in shape (321. 
If the power of the driving laser is further increased, the 
amplitude of the mechanical oscillation can be comparable 
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with the wavelength of the laser, i.e., the system reaches 
the extremely-large-amplitude regime (ELAR) J33) . In this 
regime, the optomechanical coupling should not be simply 
considered as linear or quadratic, but should be treated di¬ 
rectly as a exact function of the position of the mechanical 
resonator without approximation. In our previous work |33| , 
we studied purely classical dynamics in the ELAR of an op¬ 
tomechanical system which is a Fabry-Perot cavity with a 
perfectly reflecting and movable end mirror on one side. In 
that model, multiple optical cavity modes of different orders 
may be excited and coupled with the movable mirror via radi¬ 
ation pressure during the mechanical oscillation in the ELAR. 

In this article, we consider another optomechanical model 
in which a partially reflecting and flexible membrane is 
placed inside a Fabry-Perot cavity. Recently, this type of 
membrane-in-the-middle optomechanical system (MIMOS) 
has been widely investigated 1 34] 351, due to their superi¬ 
ority to the traditional counterparts. The optical and me¬ 
chanical functionality are segregated to physically distinct 
structures, so both excellent optical and mechanical proper¬ 
ties can be achieved simultaneously. In addition, when the 
membrane is placed at a node (or anti-node) of the intracav¬ 
ity standing wave, the system can exhibit quadratic optome¬ 
chanical coupling |!34j], this allows for direct measurement of 
the square of the membrane’s displacement, and thus quan¬ 
tum non-demolition readout of the membrane’s eigenstates. 
In this opomechanical model, the cavity resonant frequen¬ 
cies are periodic in the membrane’s position. In the ELAR, 
with an external driving laser of fixed frequency, one single 
cavity mode is excited multiple times (rather than multiple 
cavity modes are excited two times as in f33| ) and coupled 
with the membrane during one cycle of the mechanical os¬ 
cillation. In this article, we study classical dynamics of the 
system and the evolution of the quantum fluctuations around 
the classical orbits in the ELAR. We also calculate the entan¬ 
glement between the optical cavity mode and the membrane 
in the ELAR, as entanglement is an important and precious 
resource for quantum information processing tasks |36] - |43) . 

We organize the article in the following way. Sec. [^intro¬ 
duces the model and the Hamiltonian, and gives the dynami¬ 
cal equations of the system. In Sec. [3] the classical dynamics 
of the system is studied by analyzing the limit cycles of the 
mechanical oscillation of the membrane in the phase space. 
The evolution of the quantum fluctuations around the classi¬ 
cal orbits and the quantum entanglement between the optical 
cavity mode and the mechanical resonator is studied in Sec. 
[4] Finally, Sec. [5]gives a summary of this work. 


is of length L with the end mirrors fixed at x — +L/2. It is 
driven by an external laser of frequency u>i and power P. If 
the driving laser is turned off, the static equilibrium position 
of the membrane is q s . 



Figure 1 (Color online) Schematic of a MIMOS. A partially reflecting, 
flexible thin dielectric membrane is placed inside a Fabry-Perot cavity. The 
cavity of length L is driven by an external laser of frequency tu/ and power 
P. The membrane can be considered as a mechanical resonator with intrinsic 
frequency aj m , damping rate y, static equilibrium position q s and dynamic 
position q. It is coupled with a cavity mode with ^-dependent frequency 
a> c (q) and decay rate k via radiation pressure. 


In this system, radiation pressure provides the dominant 
optomechanical coupling which is typically dispersive, im¬ 
plying that the primary effect of the membrane’s motion is to 
shift the frequencies of the optical cavity modes. In the sim¬ 
plest case that \r c \ — 1 and q — 0, where r, is the reflectivity 
of the membrane and q is its position, the cavity is divided 
into two uncoupled subcavities with equal lengths L/2 by the 
membrane. The resonant frequencies of the two subcavities 
are both u>„ = hmcj L with mode number n and wavelength 
A„ = L/n. In the general case, the membrane is partially re¬ 
flecting, i.e., \r c \ 4 1, the two subcavities are coupled and 
each pair of the twofold degenerate modes of frequencies u> n 
splits into a pair of nondegenerate modes of frequencies io, ue 
and <jj n o as below ||44 j, 

c . , 

^n,e (*?) — + — sill 

TIC C 
tOn 

( 2 ) 

Thus, the frequencies of the optical cavity modes are periodic 
functions of q, as shown in Fig. If one cavity mode of 
frequency co c (q) = w„ /r (<T = e, a) and decay rate k satisfies 
the below relations: 


|r c |cos(^) 

An 


-sin ‘(|r c |), (1) 


|r c [cos(^p) 

An 


- sin '(|r f |). 


2 Hamiltonian and Dynamical Equations 


min(w„, cr ) < o)i < max(cu„,o-), 


We consider a MIMOS as shown in Fig. |T] in which a par¬ 
tially reflecting, flexible thin dielectric membrane is placed 
inside a Fabry-Perot cavity. The flexible membrane in this 
model can be considered as a mechanical resonator of intrin¬ 
sic frequency a> m , mass m and damping rate y. The cavity 


and 


nc 2c . _j .. 
k •« — - — sin (|r c |), 


where nc/L - 2c/Lsin 1 (|r c |) is the minimum frequency dif¬ 
ference between the cavity mode and the adjacent cavity 
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modes, then it is appropriate to consider only this single cav¬ 
ity mode participating the coupling with the membrane dur¬ 
ing the mechanical oscillation. 



Figure 2 (Color online) The frequencies of the cavity modes, which are 
periodic in the position of the membrane. 


In most early works, the displacement q - q s of the mem¬ 
brane is assumed to be very small. In this small-amplitude 
regime, a> c (q) can be approximately expanded to at most the 


second order 1341 around the static equilibrium position q s . 


a>c(q) w a> c0 + gi(q - q s ) + g 2 (q - q s ) 2 . 


where gi = a>' c (q s ) is the linear coupling strength and g 2 = 
u>"(q s )/2 is the quadratic coupling strength. 

If the optical cavity is driven by a high-power laser, the 
small-amplitude assumption will be no longer valid and the 
expansion of a> c (q) to the second-order will not be a good 
approximation. It is necessary to deal with oj,(q) directly 
without any approximate expansion. If the power P is high 
enough, the amplitude A of the mechanical oscillation may 
be comparable with the wavelength A; of the driving laser. 


A/Ai ~ 1, 


i.e., the system reaches the ELAR. In this case, as oj c (q) 
is a periodic function of q, each time the membrane passes 
through the positions that satisfy oj c (q) = a>i, the optical cav¬ 
ity mode is excited. So during one whole cycle of the me¬ 
chanical oscillation, the optical cavity mode is excited mul¬ 
tiple times and be coupled with the membrane via radiation 
pressure. It should be noticed that, this case is different from 
the one in our previous work. Ref. {33) , in which the cavity 
resonant frequencies are monotonic functions of the position 
of the movable end mirror, so in the ELAR, multiple optical 
cavity modes of different orders are excited twice during one 
cycle of the mechanical oscillation. Here, we assume that 
the size of the membrane is much larger than its amplitude 
of mechanical oscillation. In this situation, we can treat the 
membrane as a harmonic resonator and neglect the effect of 
mechanical nonlinearities {45). In the frame rotating at fre¬ 


quency a>i, the Hamiltonian of the system reads 


t p 2 1 2 2 

H =h[oj c (q) - u)i]a'a+ — + -mtojq - q s ) 
zm z 

+ hai(a + a' i ) + H K + H y , 


(3) 


where ci and a’ are the bosonic annihilation and creation 
operators of the optical cavity mode, q and p are the posi¬ 
tion and momentum operators of the membrane, a/, is the 
complex amplitude of the driving laser field which satisfies 
|ar/J 2 = 2 kP/Hoji. Here, without loss of generality, we set ai 
to be real. H K denotes the coupling between the optical cavity 
mode and the vacuum bath that leads to the decay rate k. H y 
refers to the interaction between the membrane and the ther¬ 
mal reservoir which is the cause of the damping rate y. The 
quantum Langevin equations (QLEs) for operators q, p and a 
can be easily derived from the Hamiltonian as follows, 


d „ p 

~vA = 
dt m 

d . 4 . 2 

—p = -hu) c (q)a'a - mojjq - q s ) - yp + f), 

d r- 

—ci = -i[u) c (q ) - a>i]d - ia l - «d + v 2 /rd,„, 
dt 


(4) 

(5) 

( 6 ) 


where f] is the mechanical Brown noise operator and d m is the 
optical vacuum input noise operator. In the presence of strong 
external driving, we can rewrite each Heisenberg operator as 
the sum of classical mean value and quantum fluctuation op¬ 
erator as below. 


q = qo(t) + q-5q, 

P = Po(t) + p-AP, 
a = a(t) + da, 

where q- = ^h/maj m and p z = ytunu> m are the zero-point 
fluctuations of the membrane’s position and momentum, re¬ 
spectively. a> c (q) and a>' c (q) in Eqs. can be approxi¬ 

mately expanded around the classical mean position quit) as 
below. 


a> c (q) = cJciqo) + oj' c (q 0 )q z 6q, 
aj’ c (q) = oj'Jqo) + aj"(qo)q?Aq, 


Thus, we can obtain the classical dynamical equations of the 
system: 


qo 


dt 
d 

dt‘ 
d 

~r a 

dt 


:Po 


Po 
m ’ 

-naj' c (q 0 )\a\ 2 - moj 2 m (q 0 - q s ) - yp 0 . 


-i[aj c (q n ) - a>i]a - ia L - ko. 


(7) 

( 8 ) 
(9) 


and the QLEs for the quantum fluctuations: 


—6q=u) m 5p, (10) 

at 

-jdp = - a> m 6q ~ q:a>”(qo)\a\ 2 dq - q z oj' c (q 0 )a*6d 
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- q z oj' c (qo)a6a^ - y6p + f, (11) 

d 

—6a = - i[oj c (qo) - u){]6d - iq : aj c (qo)adq - k5ci 

+ 'llKClin, ( 12 ) 


where f = q/q z . The noise operators and a m have zero 
mean values and are characterized by their auto correlation 
functions [461 

{a in (t)a] n (f)) = 6(t - t'), (13) 

and 


(l(0l(O) 


y_ 


f 


— p -ioKt-f) 
2 n 


U) 


coth 




where kg is the Blotzmann constant and T is the tempera¬ 
ture of the membrane. The noise operator is not delta- 
correlated and therefore does not describe a Markovian pro¬ 


cess |46|. However, quantum effects are achievable only 
when the membrane has a high mechanical quality factor 
(Q m = a> m /y » 1). In this limit, £ becomes delta-correlated 

ET 

(l(0f(O + l(f')Ko) /2 - y(2n th + l)tf(r - f), (14) 


where n t h = [exp (fia),„/kBT) - 1] 1 is the mean thermal exci¬ 
tation number, and we recover a Markovian process. 


3 classical dynamics 



Figure 3 Limit cycles in the phase space of the mechanical resonator 
scanned by qo and po with different values of P, which are: (a) P = 0.02 W, 
(b) P = 0.095 W, (c) P = 0.21 W, (d) P = 0.27 W, as specified by the vertical 
black dashed lines in Fig. [T] 


Firstly, we focus on the purely classical dynamics of the sys¬ 
tem. We integrate Eqs. ([7|-(j9]) numerically by using a fourth- 
order Runge-Kutta algorithm. Considering the experimen¬ 
tal feasibility [34 35), we set the parameters as follows, the 


intrinsic frequency, mass, damping rate and reflectivity of 
the membrane are cj m — 2n x 10 s Hz, m — 5 x 10 14 kg, 
y = 10 2 oj, h and r c = 0.8. The length of the cavity is L - 6 
cm. During the mechanical oscillation of the membrane, only 
one single cavity mode is excited, it is of order n = 60000, 
wavelength A„ = 1000 nm, frequency oj, = w„ „ and decay 
rate k = 50a> m . It should be noticed that we set the cavity 
in the unresolved sideband regime (k > u> m ) to make sure it 
can respond quickly enough to the fast mechanical oscilla¬ 
tion. The static equilibrium position of the membrane is set 
to be q s = A„/ 8. The system is driven by a external laser of 
frequency oj, = 0 J n . 0 (q s ). 



- (Is)IK 



P{ W) 


Figure 4 (Color online) (a) Definition of the average amplitude as A = 

KKn + Am ax )/ 2. (b) Attractor diagram on a plane of A and P. The vertical 
black dashed lines specify some values of P distributed in different regions. 
Limit cycles for these values of P are plotted in Fig. [3] 

From the numerical solutions of Eqs. 0-0, we plot the 
limit cycles in the phase space of the mechanical resonator 
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(C) (d) 



Figure 5 (Color online) The variation of the optical cavity occupation with the mechanical oscillation of the membrane. The values of P are: (a) P = 0.02 W, 
(b) P = 0.095 W, (c) P = 0.21 W, (d) P = 0.27 W, as specified by the vertical black dashed lines in Fig. [4] Here, we only consider the case that the membrane 
oscillating with the smallest limit cycles. 


scanned by q a and po as shown in Fig. [3] Each limit cy¬ 
cle corresponds to a stable self-sustained oscillation in the 
long-time limit. The limit cycles are in the shape of sawtooth- 
edged ellipses similar as the case in Ref. [331. However, there 
are some differences between the present limit cycles and 
the ones in Ref. |33) . In Ref. f33| , the limit cycles expand 
when the movable mirror moves forward and shrink when 
the mirror moves backward. It is because the cavity field 
is distributed on one side of the movable mirror. When the 
mirror moves forward, the radiation pressure of the excited 
modes always does positive work on it and makes it sharply 
accelerate, as a reflection in the phase space, the limit cycles 
expand. And vice versa, when the mirror moves backward, 
the limit cycles shrink. While in the present model, the cav¬ 
ity field is distributed on both sides of the membrane. When 
the membrane moves ,whether forward or backward, the ra¬ 
diation pressure of the excited cavity mode alternately does 
positive work and negative work on it. So it experiences al¬ 
ternately sharp acceleration and deceleration. As a reflection 
in the phase space, the limit cycles alternately expand and 
shrink whether along the forward direction or the backward 
direction. 


As shown in Fig. [3] there are multiple different limit cycles 
in the phase space, meaning that the mechanical resonator ex¬ 
hibits dynamical multistability. The mechanical oscillation is 
stable only when the total net work done by radiation pres¬ 


sure is balanced with the dissipative energy during one whole 
cycle. At some parameters, there are multiple stable oscil¬ 
lations that can satisfy the energy-balance condition. So the 
mechanical resonator can exhibit dynamical multistability. 

To demonstrate dynamical multistability concisely, we de¬ 
fine average amplitude as A — ^(A„ iin +Aj „ ax )/2 as shown in 
Fig. [4ja), where A min and A max are respectively the minimum 
and maximum amplitude of a limit cycle. It should be no¬ 
ticed that with the effect of radiation pressure, the dynamical 
equilibrium position qo is shifted from the static equilibrium 
position q s . However, in the EFAR, this shift is very small 
compared with the amplitude and can be neglected. So A min , 
A max , and A are all defined with the origin at the static equi¬ 
librium position q s . Fig. 0b) shows the dependence of A on 
P, which can be considered as an attractor diagram. The dis¬ 
creteness of A reveals that the energy-balance condition leads 
to an amplitude locking effect. At some parameters, the me¬ 
chanical resonator exhibits dynamical multistability, so A can 
take multiple values for a fixed P as shown in Fig. |4]b). 

We demonstrate the variation of the photon number in the 
cavity with the mechanical oscillation of the membrane in 
Fig. [5] As the frequency of the cavity mode is periodic in the 
position of the membrane, every time the membrane passes 
though the positions that satisfy u> c (qo) = op (in our paramet¬ 
ric space, that is qo = q s + kA n /4, k e N), the cavity mode 
is excited, and as a result, the photon number in the cavity 
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experiences a maximum value. 


4 quantum entanglement 


Once the classical dynamics of the mean values quit), po(t), 
and a(t) is obtained, we can solve the dynamics of the cor¬ 
responding quantum fluctuations 6q , <5/3, and 5a. To proceed, 
we introduce the position and momentum quadratures for the 
cavity mode and the corresponding input noise: 


6x = 


6y : 


6a + 6ci 1 

v^T’ 

6a - 6a 1 

iV2 ’ 

a &» + oJ„ 
V2 ’ 
a m - at. 


yin ~ 


/V2 


We write all quadratures and noise operators in terms of vec¬ 
tors: 


u - [6q,6p,6x,6y] T , 
e = [0,|, V2 Kx in , V2 Ky in ] r . 


Then, Eqs. (fl0|-([T2| can be written in compact form as 
ii = A(t)u + e, 

with the real time-dependent matrix 

0 a> m 0 0 5 


Mt) = 


-Q„,(0 -y -G x (t) -G y (t) 
Gy(t) 0 -K A (?) 

-G x (t) 0 -A(0 -k , 


(15) 


(16) 


where 

Mt) = co c (qo) - u»i 


(17) 


is the effective detuning, G x (t) and G y (t) are, respectively, real 
and imaginary parts of the effective coupling 

G(t) = yf2q z (v' c (q 0 )a = G x (t) + iG y (t), (18) 


and 

Q,n(t) = oj m + q 2 z cj"(qo)\a\ 2 . (19) 

Since H and a in are zero-mean quantum Gaussian noises and 
Eq. <[T3 is linearised, the state of the quantum fluctuations 
converges to a time-dependent Gaussian state [|48j, fully char¬ 
acterized by the covariance matrix V(t) whose elements are 
defined as 


Vtj(t) = ( Ui(t)Uj(t ) + uj(t)ui(t )) /2. (20) 


From Eq. (15 i, we can easily obtain an equation of evolution 
of V(t): 

V(t) = A(t)V(t) + V(t)A r (t) + D, (21) 


(a) 


(b) 



* - 0.15 - 0.15 

(C) ( d ) 



Figure 6 (Color online) The variation of the quantum entanglement with the mechanical oscillation of the membrane. The values of P are: (a) P = 0.02 W, 
(b) P = 0.095 W, (c) P = 0.21 W, (d) P = 0.27 W, as specified by the vertical black dashed lines in Fig. [4] Here, we only consider the case that the membrane 
oscillating with the smallest limit cycles. The temperature of the membrane is T = 1 mK. 
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Figure 7 (Color online) (a) Time evolution of qo, po, \a\ 2 and E^. (b) The respective zoom-in regions between the two vertical black dashed lines in (a). 
Here, the power of the driving laser is P = 0.21 W, the membrane oscillates with the smallest limit cycle, and the temperature of the membrane is / = 1 mK. 


where D is the diffusion matrix of the noise and is defined as: 

6(t - t')Djj(t) = ( ei (t)e/t') + /2. (22) 

From Eq.(p~3]> and Eq.([T4|, we can get that: 

D = diag[0,r(2« /A + 1 ),a-,*-]. (23) 

The quantum entanglement between the optical cavity mode 
and the membrane can be quantified by the logarithmic neg¬ 
ativity En, a quality which has been proved as a measure of 
entanglement f49) . In the case of continuous variables, E n 
can be defined as (50] 

En = max[0, - ln(2/7 )], (24) 

where 

rf = 2 _1/2 {S - [S 2 - 4det V] 1/2 } 1/2 , 

£ = det V\ + det Vi - 2 det V3, 


where V\, V 2 , and V 3 are 2x2 subblock matrices of V(t) as 



y 3 

V 2 


)■ 


By numerically integrating Eqs. (|7]i-(j9]i together with Eq. 

we can obtain the evolution of the quantum entangle¬ 
ment by using the logarithmic negativity. We show the varia¬ 
tion of the quantum entanglement with the mechanical oscil¬ 
lation of the membrane in Fig. [6] Similar as the cavity oc¬ 
cupation, the quantum entanglement also experiences a max¬ 
imum value at positions q 0 = q s + kA n /4. 


For more clearly compare the classical dynamical process 
and the evolution of the quantum entanglement, we plot the 
evolution of qo, po, |a| 2 and £\v together in the time domain, 
as shown in Fig. [7] It can be seen that, there is some synchro¬ 
nism between the classical dynamical process and the evo¬ 
lution of the quantum entanglement. Every time the mem¬ 
brane passes though the positions qo = q s + kA n /4, the fre¬ 
quency of the cavity mode a> c (qo) meets the frequency of the 
external driving lasr up. The cavity mode is excited and the 
cavity occupation experiences a maximum value. The mem¬ 
brane sharply accelerates or decelerates due to the kick ef¬ 
fect of the radiation pressure. The magnitude of the effec¬ 
tive coupling Git) = yl2q z uj' c (qo)a increases, and as a result, 
the quantum entanglement experiences a maximum value. 
While, when the membrane moves away from the positions 
qo - <ls + kA„/4, the cavity occupation decreases rapidly to 
almost zero. The membrane carries out damped harmonic 
oscillation. The magnitude of the effective coupling Git ) de¬ 
creases to almost zero, so the quantum entanglement vanishes 
due to the decoherence. 

5 Summary 

We have studied classical dynamics and quantum entangle¬ 
ment of a MIMOS in the ELAR, in which the optical cav¬ 
ity mode can be excited multiple times during one cycle of 
the mechanical oscillation of the membrane. Every time the 
optical cavity mode is excited, the radiation pressure kicks 
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the membrane and makes it sharply accelerate or decelerate. 
So the membrane can present self-sustained oscillations with 
limit cycles in the shape of sawtooth-edged ellipses. We have 
demonstrated that the membrane can exhibit dynamical mul¬ 
tistability in a wide range of parametric space. We have also 
studied the dynamics of the quantum fluctuations around the 
classical orbit and calculated the evolution of the quantum 
entanglement between the membrane and the optical cavity 
mode during the mechanical oscillation. We have shown that 
there is some synchronism between the classical dynamical 
process and the evolution of the quantum entanglement, and 
this may find some applications in quantum information pro¬ 
cessing. 


This work is supported by the National Natural Science Foundation of China 
under Grant Nos. 11175094 and 91221205 and the National Basic Research 
Program of China under Grant No. 2011CB9216002. 
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